Detecting Cosmic Strings in the CMB with the Canny Algorithm 
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Line discontinuities in cosmic microwave background anisotropy maps are a distinctive prediction 
of models with cosmic strings. These signatures are visible in anisotropy maps with good angular 
resolution and should be identifiable using edge detection algorithms. One such algorithm is the 
Canny algorithm. We study the potential of this algorithm to pick out the line discontinuities 
generated by cosmic strings. By applying the algorithm to small-scale microwave anisotropy maps 
generated from theoretical models with and without cosmic strings, we find that, given an angular 
resolution of several minutes of arc, cosmic strings can be detected down to a limit of the mass per 
unit length of the string which is one order of magnitude lower than the current upper bounds. 

PACS numbers: 98.80.Cq 



INTRODUCTION 



Cosmic strings [l[ are one-dimensional topological de- 
fects which arise during phase transitions in the very 
early universe. Since they carry energy, they will lead to 
density fluctuations and cosmic microwave background 
(CMB) anisotropies (see e.g. [|J for reviews on cosmic 
strings and structure formation). Causality implies that 
the network of strings which forms during the phase tran- 
sition contains infinite strings. Once formed in the early 
universe, the network of strings will approach a "scal- 
ing solution" which is characterized by of the order one 
infinite string segment in each Hubble volume, and a dis- 
tribution of cosmic string loops which are the remnants 
of the previous evolution. In particular, this implies that 
in any theory which admits cosmic strings, a network of 
strings will be present during the time period relevant 
to the CMB, namely between the time of recombination 
t rec and the present time irj. 

A network of cosmic strings will generate a scale- 
invariant spectrum of cosmological perturbations [3| . As 
first discussed by Kaiser and Stebbins (KS) [H, the non- 
Gaussian nature of the density field produced by strings 
will lead to a distinctive signature (which we will call KS 
signature) in CMB anisotropy maps, namely line discon- 
tinuities. These line discontinuities are a consequence of 
the non-trivial nature of the metric produced by a cosmic 
string: space perpendicular to a cosmic string is a cone 
with deficit angle given by 



a = 8nGfi . 



(1) 



where \i is the mass per unit length of the string and 
G is Newton's gravitational constant Q. Since cosmic 
strings have a tension comparable to their mass per unit 
length, they will typically be moving with a relativistic 
transverse velocity v. As illustrated in Figure 1, if we are 
looking at the CMB in direction of the string, we will see 
the photons passing on different sides of the string with 
a Doppler shift 

— = 8ir-f(v)vGn, (2) 
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FIG. 1: Geometry of the Kaiser-Stebbins effect: Photons 
passing on the two sides of the moving cosmic string obtain 
a relative Doppler shift for the observer who is at rest. 



where 7(f) is the relativistic gamma factor. Looking in 
direction of the string, we will see a line in the sky across 
which the CMB temperature jumps by the above amount. 

Any cosmic string which photons reaching us today 
pass on their way from the last scattering surface will 
produce an effect. The most numerous strings, however, 
are those present close to the time of last scattering t rec . 
At that time, the typical curvature radius of a long string 
segment is of the order of the Hubble radius at t rec . The 
corresponding angular scale today is about 1°. In order 
to be able to identify the cosmic signature as a line dis- 
continuity, as good an angular resolution as possible is 
required. If probed with an angular resolution not signif- 
icantly smaller than 1°, the KS signature will be washed 
out. In this work we will study the potential of surveys 
with an angular of several minutes of arc to detect the 
KS signature for strings. Less important than the angu- 
lar resolution of the survey is the total sky coverage, as 
long as it includes a sufficiently large number of Hubble 
patches at t rec . Obviously, increasing the sky coverage 
will reduce the standard deviation of the results and thus 
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lead to some better discriminating power. 

There has been surprisingly little work devoted to de- 
tecting the KS signature. An early study @ showed that 
the angular resolution of the WMAP satellite would not 
be adequate to pick out the KS signature, even for val- 
ues of G/i for which the cosmic strings would dominate 
the power spectrum of density fluctuations. In order to 
detect the KS signature, new statistical methods must 
be employed. In 0, a matched filtering method was 
applied to the WMAP data to search for cosmic string 
signatures. In [8], new statistics (such as a statistic 
measuring the connectedness of neighboring temperature 
steps or a decomposition of the temperature map into 
constant background temperature, Gaussian noise plus 
straight string step discontinuities) were introduced and 
applied to the WMAP data to search for strings. Based 
on the null results of both studies, an upper bound on 
G/i of G/i < 10~ 6 could be set. However, this bound 
is not competitive with other existing bounds on the 
string tension coming from the observed acoustic peak 
structure of the CMB angular power spectrum [9j (which 
yields G/i < 1CP 7 ) and from pulsar timing measurements 
[l0| (which yield limits of between G/i < 2 x 10~ 8 and 
G/i < 10~ 5 ' 5 , the difference being due to conflicting ways 
[ll| in which the pulsar timing data is analyzed to ob- 
tain limits on the amplitude of the stochastic background 
of gravitational waves, and due to different assumptions 
about the distribution of cosmic string loops, the main 
source of gravitational radiation from a string network 
- see [HJ for a recent discussion). Note that an advan- 
tage of using the KS signature over pulsar bounds is that 
the KS signature is more robust (less dependent on the 
unknown distribution of cosmic string loops) than the 
pulsar bounds (which depend sensitively on the details 
of the cosmic string loop distribution) [32| 

As discussed above, we expect to be able to set much 
stronger limits on G/i from searches for the KS signature 
by looking at small-scale anisotropy maps, provided we 
find a good statistic to identify line discontinuities. The 
goal of this paper is to explore the potential of one spe- 
cific edge detection technique, the Canny algorithm [13[ 
to pick out these line discontinuities in CMB anisotropy 
maps. 

One reason for the relative lack of previous work on de- 
tecting the KS signature of cosmic strings is the fact that 
interest in cosmic strings as a source for structure in the 
universe [33[ decreased dramatically after the discovery 
of the sharpness of the acoustic peaks in the CMB angu- 
lar power spectrum [lH [l5T | . Since perturbations seeded 
by strings are not coherent [l6[ , a scenario with only cos- 
mic strings as the source of density perturbations would 
only produce one broad Doppler peak rather than the 
observed narrow acoustic oscillations (l7l [H, [l9| . 

However, recent years have seen a resurgence of interest 
in cosmic strings. This was fueled by two developments. 
Firstly, it was realized [201 ] that in many supersymmetric 
particle physics models, cosmic strings are formed after 
inflation, and thus contribute to but do not completely 



replace inflationary perturbations as the seeds for struc- 
ture formation. Secondly, it has recently been realized 
that models with cosmic superstrings [211 ] may well be 
viable [22| . They could, for example, be generated as the 
remnant of brane annihilation processes in brane infla- 
tion models [23j ], or they may play an important roles 
in inflationary models in warped backgrounds [24] . Cos- 
mic superstrings may also be left behind after the initial 
Hagedorn phase in string gas cosmology [25| , where they 
would add an additional component to the spectrum of 
fluctuations produced by thermal string gas fluctuations 
[2^ ]. Thus, it is of great interest to find new ways to 
search for signatures of strings in cosmological data. 

The outline of this paper is as follows: In the follow- 
ing section we review the Canny algorithm and describe 
its application to our problem. Section 3 discusses how 
the temperature maps with and without strings are pro- 
duced, introduces the parameters chosen in the specific 
simulations, and presents our results. We conclude with 
a summary and a discussion of the caveats of the current 
analysis and prospects for future work. 



II. THE CANNY ALGORITHM 

The Canny algorithm was developed in 1986 as a tech- 
nique to detect edges in images [13| such as the two- 
dimensional images considered here. When applied to a 
map of raw data, the algorithm is intended to produce a 
map tracing the edges in the map, the lines across which 
the intensity contrast is largest. 

The first step in the algorithm is to filter the data to 
eliminate point source noise. The filtering is achieved us- 
ing a convolution of the map with a Gaussian filter. The 
filtering length is a free parameter in the algorithm. It 
must be chosen sufficiently large to eliminate unwanted 
point source noise, but must be smaller than the char- 
acteristic size of the structure in the maps which one is 
trying to identify. Let us denote the original map by 
M(i,j), and the filter by F(i,j). The filtered map FM 
is then given by 

FM(i,j) = ^M(i-k,j-l)F(k,l), (3) 
k,i 

where the pixel points of the map are denoted by the 
labels In our simulations, the filtering length is 

taken to be 1.5 grid units. 

The second step of the algorithm is to find the gradient 
of the filtered image, a vector obtained by taking the dis- 
crete derivative of the map in each of the two directions. 
The gradient vector at a grid point (i,j) is denoted by 
(G x , G y ), where G x and G y are the derivatives in the re- 
spective directions. The edge strength \G\ at any given 
grid point is defined by 

\G\ = yJ\G x \* + \G y \*, (4) 
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and the edge direction is given by the angle 

9 = arctan{G y /G x ) . (5) 

Since we are working on a grid, it makes sense to replace 
the edge direction by one of the eight distinguished di- 
rections on a quadratic grid, namely the four directions 
along the coordinate axes and the four diagonal direc- 
tions. 

The average maximal gradient in maps without cosmic 
strings (the average taken over all of our runs) is denoted 
by G m . For a specific run, the maximal gradient is ob- 
tained by scanning the entire map and searching for the 
maximal value of the grid strength. The value of G m is 
used to set the thresholds discussed below. 

The next step of the Canny algorithm is to find grid 
points with gradients which are maximal when we vary 
the point in direction of the gradient. Grid points which 
are not local maxima of the gradient are assigned a num- 
ber 0. 

At this stage, two thresholds must be set, an upper and 
a lower threshold t u and ti . The algorithm first finds local 
maxima in the gradient along one of the four directional 
axes on the grid. Next, it goes through each of the list 
of local maxima and determines whether the gradient is 
greater than the upper threshold fraction of G m , i.e. 

\G\ > t u G m . (6) 

Grid point thus selected are marked with a number 1 . If 
the gradient is larger than the lower threshold fraction of 

Gjn 5 

\G\ > UG m , (7) 

the grid point is assigned a value of 1/2. For points 
marked with 1 or 1/2, the direction of the gradient is 
also stored. For each grid point with a value of 1/2, the 
algorithm next checks whether there is a grid point with 
value of 1 or 1/2 which is in a direction perpendicular 
to the gradient and whose gradient is parallel or next to 
parallel (thus three possible directions in total) to the 
initial gradient. If such a grid point is found, the algo- 
rithm continues until it either does not find another point 
satisfying the criteria or else finds a point satisfying the 
criteria marked with a 1. If it finds such a point, then all 
of the points (initially marked 1/2) found along the way 
are marked as 1. Neighboring (in all eight directions) 
grid points marked by 1 are then said to belong to the 
same edge. These are edges across which the gradient of 
the map is large. 

The reason for using two thresholds is as follows: we 
want the algorithm to find lines with a large gradient. 
At the same time, we do not want noise to lead to points 
along a line of large gradients to be missed (and thus the 
lines cut) just because noise has reduced the amplitude 
of the edge strength at one point along the line. 

In this way, the algorithm produces a list of grid points 
belonging to edges, and this list can be read out as a 



map of edges. In our simulations, the upper and lower 
thresholds are t u — 0.5 and ti — 0.4. 

It is important to quantify the edge maps. A simple 
way to do this is to produce a histogram of edge lengths. 
Thus, the implementation of the Canny algorithm is ar- 
ranged such as to output a list of edge lengths which can 
then be used to produce a histogram of edge lengths. The 
histogram contains useful information about the presence 
of edges in the map. 

III. APPLYING THE CANNY ALGORITHM TO 
CMB ANISOTROPY MAPS 

Eventually, we would like to apply the Canny algo- 
rithm to actual CMB maps. For this initial feasibility 
study, however, we will apply the algorithm to theoreti- 
cal maps produced in numerical simulations. We wish to 
compare temperature maps for the "standard" Gaussian 
ACDM model with maps in which a network of strings 
contributes a fraction / of the total power. The maps are 
characterized by the angular scale of the survey (we take 
the survey area to be square) and by the angular scale of 
the grid, i.e. the angular resolution of the survey. 

The Gaussian maps are produced in the following way: 
We start with the angular power spectrum C; of cosmic 
microwave anisotropics taken from the CMBFAST sim- 
ulation (27| with the appropriate cosmological parame- 
ters [34] ]. Since we have in mind applications of our 
algorithm to small angular scale surveys, we perform 
a "flat sky" approximation [28|. We introduce a two- 
dimensional Cartesian coordinate system covering the 
survey area which we take to be rectangular, choosing 
the upper right corner of the survey area to be the origin 
of the coordinates. We need to compute the tempera- 
ture field Tq(x,d) of the map at the coordinate values 
(x, y) corresponding to the grid points. This map is de- 
termined by an inverse fast Fourier transform from the 
temperature values T(k) in Fourier space. 

For each vector k in Fourier space, we find the integers 
l\(k) and h(k) which bracket the I value l{k) correspond- 
ing to k and take a linear interpolation of the values of 
Ci . Let us denote the result of this linear interpolation 
by CW), where k = \k\. The value of T(k) is then given 
as follows: 

T(k) = g{k)^C^, (8) 

where g is a random variable drawn from a probability 
distribution with variance 1. 

In the above analysis, we are singling out one grid 
point to be special, namely the origin. This introduces 
unwanted phase correlations which we overcome by su- 
perimposing the results from four separate simulations, 
indicated by Ti, i = 1, ..,4: 

Ta{x,y) = ~(Ti(x,y)+T 2 (x m -x ) y m -y) (9) 
+T 3 (x m -x,y) + T 4 (x, y m - y)) , 
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where x m and y m are the maximal x and y values of the 
survey volume (the pre- factor of 1/2 is required in order 
to maintain the original standard deviation). 

A network of cosmic strings will produce a temper- 
ature map Tcs{x,y). In the following, we take a toy 
model for a temperature map produced by strings intro- 
duced by Perivolaropoulos [a [2!| • The toy model takes 
into account that at all times between t rec (the time of 
recombination) and the present time £o ? the network of 
strings is described by a scaling solution in which there 
are a fixed number of long string segments (strings which 
are not loops with radius smaller than the Hubble radius) 
crossing each Hubble volume. Each such string gives rise 
to a Kaiser-Stebbins line discontinuity ([2]) in the temper- 
ature map. The cosmic string network is continuously 
evolving via motion of the strings and string interactions 
which lead to the production of string loops. Thus, in 
each Hubble time step the string network can be taken 
to be uncorrelated. 

In the algorithm, we divide the time interval t rec < 
t < to into 15 Hubble time steps. In each time step t we 
lay down a network of strings at random, uncorrelated 
with the network at the previous time step. We take the 
network to consist of straight string segments of length 
ait. In order to avoid missing strings at the edges of the 
survey area, the survey area is extended in each direction 
by a Hubble distance. The program runs through all 
points of the simulation volume and picks points to be 
midpoints of a string segment with a probability chosen 
such that the average number of strings in the Hubble 
volume equals the number N of the scaling solution. The 
directions of the strings are chosen at random, as are the 
string transverse velocities. To take into account also 
the projection onto the last scattering surface, we add a 
temperature 

St i 

— = -irSirGv, (10) 

to one side of the string projection onto the last scattering 
surface, and subtract it from the other (so as to maintain 
the average temperature of the CMB). Here, v stands for 
the maximal value of v~/(v), where v is the transverse 
velocity of the string, and "f(v) is the relativistic 7 factor 
associated with v. Also, r is a random number between 
and 1, to take into account both the distribution of 
velocities, and also projection effects (the formula ^fy for 
the line discontinuity in temperature is modulated by an 
angular factor if the string velocity is not perpendicular 
to our line of sight to the string) . The regions affected by 
the temperature fluctuation are rectangles on cither side 
of the string. The depth of these rectangles in direction 
transverse to the direction of the string is taken to be a 
fraction / of the Hubble radius. The length of the string 
segments is taken to be a,\ times the Hubble radius. The 
survey area is always taken to correspond to the central 
region of the simulation area. 

The free parameters of the string simulation are G/x, 
the number of strings N per Hubble volume, the length 



coefficients ot\ and /, as well as v. In our work, we will 
fix N = 10, ai = 2, / = 1 and v = 0.15. These values 
are representative of the results of numerical simulations 
of cosmic string evolution (see Q for reviews). However, 
it should be kept in mind that the results for these pa- 
rameters obtained using different numerical codes differ 
significantly. 

Let us make a couple of comments on the expected de- 
pendence of the results on the value of N. As N initially 
increases from zero, the signature of the strings should 
increase since there are more strings. However, once N 
rises above a certain critical value, then further increas- 
ing the number N will render the string distribution more 
Gaussian by the Central Limit Theorem and will hence 
decrease the sensitivity of the Canny algorithm. We ex- 
pect that the critical value of iV" will be the one for which 
the regions of the sky effected by individual string seg- 
ments begin to overlap. If the strings have length com- 
parable to the Hubble radius, the the critical value of N 
is expected to be less than N = 10. 

Increasing the string velocity v will boost the line dis- 
continuites in the temerature maps and will make our 
algorithm more effective. Since the weight in our statis- 
tical analysis of the edge maps comes from fairly short 
segments, our results are not very sensitive to the specific 
values of / and a%. 

We are interested in using the Canny algorithm to find 
signatures for cosmic strings in maps which contain both 
Gaussian fluctuations from inflation and a cosmic string 
component (which is sub-dominant in terms of its contri- 
bution to the power spectrum). The total map Tx{x,y) 
is given by the superposition of a pure Gaussian noise 
map and the map produced by a distribution of strings: 

Tr(x,y) = T G (x,y) +T C s(x,y) , (11) 

and where the amplitude of the Gaussian term is adjusted 
such that the total angular power spectrum agrees with 
the COBE results [H- 

The output of the Canny algorithm is a map of edges 
which have been picked out. In the presence of cosmic 
strings, we would expect more longer edges than in the 
absence of strings. In order to test for this effect, we 
produced histograms of the distribution of edge lengths 
for each simulation. Since both the Gaussian maps and 
the maps with cosmic strings are produced by a Gaussian 
random process (the phases of the Fourier modes for the 
Gaussian maps are picked at random, and for the string 
maps the locations of the string centers, their directions 
and their transverse velocity vectors are all random), we 
ran many (N = 50) simulations for both the Gaussian 
map and the maps with strings. This gave us statistical 
error bars on the histograms of edge lengths. In turn, this 
allows us to assign statistical weight to the difference in 
the histograms. Our results are based on this analysis. 

In our simulations, we consider the "Gaussian noise" 
maps to be given by a standard ACDM model with adi- 
abatic fluctuations with a spectral index of n s = 0.99, 
and with parameters Qb = 0.046 (baryon fraction), 



FIG. 2: Map of the CMB temperature in a 15 2 square of the 
sky for a ACDM simulation with the parameters described in 
the text. 



FIG. 3: Corresponding map of the CMB temperature in a 
simulation which includes cosmic strings with a mass per unit 
length parameter given by Gfi = 1CP 7 (and other parameters 
as described in the text. 



Qcdm — 0.224 (cold dark matter fraction), £1a = 0.730 
(cosmological constant contribution) , and Ho = 72 (Hub- 
ble expansion rate). No massless neutrinos, standard re- 
combination history and a Helium fraction of Y = 0.24 
were assumed (in the C\ spectra which were used to con- 
struct the maps). 

The parameters of the string simulation are indicated 
above. Our maps have angular extent 15 degrees by 15 
degrees. 

We fix the angular resolution and search for the min- 
imal value of the cosmic string mass parameter G/Lt for 
which the maps with the strings can be distinguished 
from the pure Gaussian maps at a statistically signifi- 
cant level. We demand significance at the 3tr level. 

Figure 2 shows a CMB temperature fluctuation map 
of a simulation without strings, and Figure 3 is the cor- 
responding map of a simulation which includes strings 
according to the prescription described above. 

The output of the Canny algorithm applied to the 
above maps is shown in Figures 4 and 5, respectively. 

Based on these results, our algorithm evaluates the dis- 
tribution of edge lengths and checks whether the distri- 
butions of the histograms with and without strings are 
statistically different. Maps with cosmic strings showed, 
as expected, a slight excess of edges of all lengths. This 
excess is more pronounced for edge lengths which lie in 
the non-Gaussian tail. 

The results for an angular resolution of 8 are presented 
in Table 1. In this table, the rows indicate the length L of 
the edge. The entries in the columns are the mean num- 




FIG. 4: Output map of the Canny algorithm showing the 
edges in the ACDM map of Figure 2. 
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FIG. 5: Corresponding edge distribution in the map of Figure 
3. 

ber [standard deviation of the mean] of edges of the maps 
which have the respective length. The first column of 
numbers is for simulations without strings, the next two 
are for simulations including a network of strings with 
the indicated mass per unit length parameter. Each his- 
togram was produced from 50 runs. From the histogram 
it follows that the difference between the distributions is 
significant (at the 3cr level) if G[i = 3.5 x 10~ 8 , but not 
if GfJL = 3 x 1CT 8 . 

For a given angular resolution, we varied the string 
mass parameter G/j, to find the limiting value below which 
the differences in the histograms ceased to be statis- 
tically significant (using Fisher's combined probability 
test). Our results are summarized in Table 2. 

Based on the above results, prospects for being able 
to use the Canny algorithm to significantly improve the 
limits on the cosmic string mass parameter G\x (or de- 
tect cosmic strings) are excellent. The angular resolution 
of the Planck satellite experiment will be about 5 , the 
South Pole Telescope is aiming for a resolution of 1 with 
a field of view of 4, 000 square degrees, and the resolution 
of the ACT telescope in Chile will be 1.5 . 



IV. CONCLUSIONS 

We have suggested a new way of looking for the spe- 
cific signature of cosmic strings, namely the Kaiser- 
Stebbins line discontinuities, in small-scale cosmic mi- 
crowave anisotropy maps. Our method makes use of the 
Canny algorithm, an edge-detection technique used often 



TABLE I: Histograms of the distribution of edge lengths L. 







no strings 


Gfi = 3 x 10~ B 


Gfi = 3.5 x 10~ 8 


L 


= 2 


165.0 [2.84] 


169.9 [3.60] 


178.0 [3.84] 


L 


= 3 


40.0 [0.86] 


42.3 [1.19] 


43.6 [1.18] 


L 


= 4 


11.8 [0.51] 


13.4 [0.58] 


14.0 [0.57] 


L 


= 5 


4.5 [0.31] 


4.8 [0.32] 


5.5 [0.34] 


L 


= 6 


1.7 [0.23] 


2.1 [0.21] 


2.2 [0.21] 


L 


= 7 


0.66 [0.13] 


0.74 [0.11] 


0.84 [0.12] 


L 


= 8 


0.18 [0.07] 


0.38 [0.09] 


0.28 [0.07] 


L 


= 9 


0.12 [0.06] 


0.1 [0.04] 


0.08 [0.05] 


L 


= 10 


0.04 [0.03] 


0.08 [0.04] 


0.04 [0.03] 



TABLE II: Critical values of the string mass parameter G\i. 



angular resolution 


Gfi 


10 


4.3 x 10" 8 


9 


4.0 x 10" 8 


8' 


3.1 x 10" 8 



in pattern recognition. 

Since the cosmic microwave temperature anisotropy 
maps induced by a network of cosmic strings are domi- 
nated by the strings present at the time of last scattering 
(when the Hubble radius, which is comparable to the cor- 
relation length of the string network at that time, is of 
the order of one degree), good small-scale angular reso- 
lution is essential in order to be able to detect strings. 
An angular resolution substantially less than 1° is re- 
quired, otherwise the signals from the line discontinuities 
are washed out. 

We have constructed CMB anisotropy maps which 
correspond to having both Gaussian "noise" with a 
nearly scale-invariant power spectrum from inflation and 
anisotropies produced by a distribution of straight string 
segments. Based on our numerical simulations, we find 
that for an angular resolution of CMB maps of 8 the 
Canny algorithm has the potential to detect strings with 
a mass per unit length /j, above a value of Gfi ~ 4 x 10~ 8 , 
close to an order of magnitude better than current limits 
based on the CMB. 

A drawback of our work is that it is based on toy model 
cosmic string CMB maps which are obtained by super- 
imposing idealized line discontinuities of straight string 
segments, not from actual string networks. Actual string 
networks contain both infinite strings and string loops. 
The infinite strings are not completely straight, but have 
a curvature radius comparable to the Hubble radius. An 
improved analysis should start from a numerical simula- 
tion of the distribution of cosmic strings, calculate the 
induced temperature anisotropies taking into account of 
all effects following the formalism set out in [19j], and 
then apply the Canny algorithm to the resulting maps. 
After completion of this work, a paper appeared [3l[ in 



7 



which CMB temperature maps on scales similar to the 
ones we are considering were constructed based on full 
cosmic string simulations. This work confirmed that the 
Kaiser-Stebbins from long string segments is the domi- 
nating visible effect, and suggested, in agreement with 
the message of our work, that strings should be visible in 
high resolution small-scale CMB anisotropy maps. How- 
ever, no string-specific statistical analyses of the maps 
like the one we are proposing were performed in [3lj . 

Another important issue which remains to be analyzed 
is the effect of instrumental noise. We have done prelimi- 
nary work on this topic and modelled instrumental noise 
by a component to the Ci spectrum which rises rapidly 
as a function of I at the scale of the angular resolution of 



the survey [351 ] . becoming dominant at the angular reso- 
lution scale. Initial results show that the efficiency of the 
Canny algorithm is not reduced. We plan to study this 
question in more detail. 
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